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CS| , Abstract The RESIK instrument on the CORONAS-F spacecraft obtamed solar 

^ ■ flare and active region X-ray spectra in four channels covering the wavelength 

O I range 3.8-6.1 A in its operational period between 2001 and 2003. Several highly 

ionized silicon lines were observed within the range of the long-wavelength chan- 
nel (5.00-6.05 A). The fluxes of the Si xiv Ly-/3 line (5.217 A) and the Si xiii 
Is^ — ls3p line (5.688 A) during 21 flares with optimized pulse-height analyzer 
settings on RESIK have been analyzed to obtain the silicon abundance relative 
to hydrogen in flare plasmas. As in previous work, the emitting plasma for each 
spectrum is assumed to be characterized by a single temperature and emission 
measure given by the ratio of emission in the two channels of GOES. The silicon 
^ ■ abundance is determined to be ^(Si) = 7.93±.21 (Si xiv) and 7.89±.13 (Si xiii) 

Oh, on a logarithmic scale with H = 12. These values, which vary by only very small 

amounts from flare to flare and times within flares, are 2.6 ± 1.3 and 2.4 ± 0.7 
^ , times the photospheric abundance, and are about a factor of three higher than 

t/3 ' RESIK measurements during a period of very low activity. There is a suggestion 

^ that the Si/S abundance ratio increases from active regions to flares. 



Keywords: Flares, Spectrum; Spectral Line, Intensity and Diagnostics; Spee- 
^ \ trum. X-ray 

0^ , 1- Introduction 

(N 

I REntgenovsky Spektrometr s Izognutymi Kristalami (RESIK: [Sylwester et al.\ 

(|2005p ) was a solar X-ray crystal spectrometer onboard the CORONAS-F space- 
' craft operating from the time of launch (31 July 2001) until a spacecraft failure 

terminated the mission in May 2003. The instrument and calibration are de- 
scribed by [Sylwester et aL|([M)5) . RESIK's wavelength range, 3.40-6.05 A, in- 
, eluded spectral lines emitted by highly ionized elements in flare and active region 
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spectra, and also continuum, which unhke many previous crystal spectrometers 
can be separated from a fluorescence background. A relatively precise absolute 
calibration makes RESIK ideal for the determination of element abundances with 
appropriate assumptions for the temperature structure of the emitting plasma. In 
this article, our motivation has been the determination of the silicon abundance 
for flares observed by RESIK, and comparison with previous estimates of coronal 
(both flare and active region) and photospheric silicon abundances. 

In previous articles reporting the analysis of RESIK spectra (see |Sylwester e^"am 
120121 and references cited therein), we have estimated the abundances of K, 
Ar, CI, and S from emission lines of helium-like or hydrogen-like ions of these 
elements during flares, using a method that involves assigning a "characteristic 
temperature" for the emitting plasma equal to that given by the ratio of the two 
channels of GOES. This method is applied here for the Si xiv Ly-/3, and Si xiii 
w3 lines. Although there are other lines in this channel due to Si xiii or Si xiv 
transitions, we defer the analysis of their fluxes to later work since the atomic 
data included in CHIANTI are at present under discussion. In Section [5] a brief 
description of the instrument is given, and in Section [3] the RESIK spectra and 
their analysis are discussed. In Section |4] the silicon abundance is derived using 
methods used in our previous work. Comparison with other silicon abundance 
estimates is made, and the atomic data used for the Si X-ray lines are critically 
examined, in Section [5] 



2. RESIK Instrument 

RESIK consisted of two spectrometers, each containing a pair of crystal-detector 
combinations to give a total of four channels to cover the wavelength range. The 
crystals were silicon or quartz, having low atomic numbers so that a background, 
formed by fluorescence of the crystal material by solar X-rays, was relatively 
small. This is an important advantage over previous spectrometers such as the 
Bragg Crystal Spectrometer on Yohkoh, operational in the 1990s, which had high- 
Z germanium crystals. The wavelength ranges of RESIK's four channels were 
3.40-3.80 A (Channel 1), 3.83-4.27 A (2), 4.35-4.86 A (3), and 5.00-6.05 A 
(4), with a pair of curved Si (111) crystals for Channels 1 and 2 and a pair of 
quartz (1010) crystals for Channels 3 and 4. These ranges included emission lines 
of hydrogen-like and helium-like Si, S, CI, Ar, and K and several groups of dielec- 
tronic satellite lines. The instrument was uncoUimated to maximize its sensitivity 
(the quoted wavelength ranges are for on-axis sources in the first diffraction 



order). A careful assessment of instrument parameters (Sylwester et ai, 2005) 
has enabled effective areas to be estimated with a precision of ~ ±20 %. Pulse- 
height analyzers on RESIK allowed the separation, or partial separation, of solar 
photons from the crystal fluorescence background. Adjustments of pulse-height 
analyzers over the course of the mission resulted in optimum settings in each 
channel by 24 December 2002, after which the fluorescence background was 
entirely removed from Channels 1 and 2 and much reduced for Channels 3 and 
4. 
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3. RESIK Observations and Analysis 



The RESIK spectra analyzed here were taken during the course of 21 flares 
that, with one exception, were observed after the optimization of the pulse 
height analyzers on 24 December 2002 so that the fluorescence background in 
Channels 3 and 4 was reduced as much as possible. The exception is a flare that 
was observed on 12 November 2002 when the pulse-height analyzers were by 
chance set to their optimum values for Channels 3 and 4. A total of 1822 spectra 
were recorded, with data-gathering intervals ranging from two seconds to a few 
minutes. The data-gathering intervals were inversely related to the incident X- 
ray fluxes to avoid detector saturation or very low photon count rates with poor 
statistical significance. Tablc[T]gives details of the flares with numbers of spectra 
collected in each period of observation. For each of these spectra, temperature 
and emission measure from the emission ratio of the two GOES channels, Tqoes 
and EMqoes (equal to N^V with = electron density and V the emitting 
volume), were evaluated. Figure [1] (upper panel) shows all 1822 Channel 3 and 4 
spectra on a colour scale stacked in order of IgoeSj with the temperature scale 
shown on the left, while Figure [T] (lower panel) shows five sample Channel 3 
and 4 spectra with the values of Tqoes indicated. This figure illustrates how 
the relative fluxes of the principal lines vary with temperature in accordance 
with their G(T) or contribution functions, equal to emission per unit emission 
measure. Thus, in Channel 4 spectra, the S XV lines near 5 A (temperature 
of maximum G(T), T^ax « 16 MK) and Si xiv Ly-/3 line (r,nax « 16 MK) 
are prominent at relatively high temperatures but the Si xiii w3 line emission 
(^max ~ 10.5 MK) extends to slightly lower temperatures. 

RESIK observations of the Si xiv Ly-/3 line at 5.217 A were analyzed by 
estimating for each spectrum the flux in the wavelength interval 5.185-5.241 A 
and subtracting a neighbouring portion of background, which includes both 
solar continuum radiation and crystal fluorescence. The values of Tqoes a-nd 
EMqoes at the times of each spectrum were derived from standard software in 
the SolarSoft IDL package, and the line flux divided by EMqoes plotted against 
Tgoes; see the left panel of Figured] The points at temperatures above « 7 MK 
have high statistical quality, but the weakness of the line at lower temperatures 
prevents measurements at Tgoes ^ 5 MK. The continuous curves in the figure 
are the theoretical contribution functions G(Te) from the CHIANTI code, defined 

by 



_ jV(Si+") iV(Si+") 7V(Si) iV(H) 
^ "^~7V(Si+") iV(Si) 7V(H) ' ^ ' 

where A''(Si^"') is the population of the excited (zth) level of the ion Si^" (n = 12 
or 13), A''(Si^")/iV(Si) is the ion fraction as a function of electron temperature 
[Te] , taken in our case from |Bryans, Landi, and Savin] (|2009p . Other quantities in 
Equation (1) are AiQ (transition probability from level i to the ground state), and 
N{}l)/Ne, equal to 0.83. For the abundance of Si relative to H, [iV(Si)/iV(H)], 
we chose as comparison to the RESIK observed points both the photospheric 
value given by ]Asplund et aL](]^DDg]) . viz. N(Si)/N(H) = 3.2 x 10"^ (equivalent 



SOLA: RESIK_Si_abund.tex; 13 December 2012; 2:06; p. 3 



B. Sylwcster et al. 



Table 1. Flares Analyzed for Si Abundance. Note: The flare indicated by " is the first of 
two overlapping flares. 



No. Date (2002/2003) Time [UT] of GOES Time Range of Number of 
Flare Peak Importance Spectra Spectra 



1 


12 Nov. 


07:50 


C5.4 


07:45 - 08:31 


93 


2 


25 Dec. 


12:07 


C3.5 


lz:Uo — 12:oo 


64 


3 


26 Dec. 


08:35 


CI. 9 


08:26 - 08:51 


53 


4 


7-8 Jan. 


23:30 


M4.9 


23:13 - 01:07 


113 


5 


09 Jan. 


01:39 


C9.8 


r\r\ 1 A r\o r\o 

00:14-02:08 


297 


6 


21 Jan. 


02:28 


08.1 


0z:2i — Oz:oo 


by 


7 


21 Jan. 


02:50 


C4.0 


02:35 — 02:44 


25 


8 


21 Jan. 


15:26 


Ml. 9 


11:57-18:56 


290 


9 


25 Jan. 


18:55 


C4.4 


18:35-19:10 


52 


10 


01 Feb. 


09:05 


Ml. 2 


08:22-10:14 


133 


11 


06 Feb. 


02:11 


C3.4 


02:01-02:35 


34 


12 


14 Feb. 


02:12 


C5.4 


01:46-02:48 


44 


13 


14 Feb. 


05:26 


C5.6 


05:00-05:57 


64 


14 


14 Feb. 


08:47 


C1.4 


08:42-09:02 


23 


15 


14 Feb. 


09:18 


05.9" 


09:11-09:57 


42 


16 


15 Feb. 


08:10 


C4.5 


06:31-11:30 


299 


17 


21 Feb. 


19:50 


C4.3 


19:46-20:10 


32 


18 


22 Feb. 


04:50 


B9.6 


04:40-05:35 


19 


19 


22 Feb. 


09:29 


C5.8 


08:29-10:15 


29 


20 


11 Mar. 


05:50 


B7.3 


04:22-07:23 


15 


21 


17 Mar. 


19:05 


XI. 5 


18:53-19:45 


32 



to A{Si) = 7.51) and the coronal value given by |Feldnian and Lamingl (j2000p . 
iV(Si)/iV(H) = 1.26 X 10-4 (A(Si) = 8.10). The red continuous curve in Figured] 
is based on the coronal Si abundance, the blue dashed curve on the photospheric 
Si abundance. Figure [2] shows that there is a close correspondence in the shape 
of the theoretical G(Te) curves with that defined by the RESIK points, but with 
a slight tendency for the observed points to be higher than the theoretical curves 
for Tgoes ^ 7 MK. This discrepancy is discussed further in Section [SJ 

The Si XIII wi line at 5.688 A is the most intense line in RESIK's Channel 4 
over temperatures ranging from non-flaring active regions ( [Sylwester, Sylwester, and Phillips, 2010l )[ 
to flares with Tgoes up to « 10 MK. The line flux, taken to be the emission 
in a band 5.66-5.71 A with a portion of the background from a neighbouring 
band subtracted, is easily measured over this range though at higher temper- 
atures there is extra, unexplained emission between the w3 line and a line 
feature peaking at 5.82 A due to several Si xii dielectronic satellites (transitions 
ls^2p^P - Is2p3p ^D: [Phillips et aJ] ^UU^). Figure |3] (left panel) shows the 
Si XIII w3 line emission divided by EMqoes against Igoes, with theoretical 
G(Te) curves shown, as in Figure[21 for photospheric and coronal Si abundances. 
The RESIK points follow the same dependence on temperature as the theory 
curves for Tgoes up to w 11 MK but attain a maximum that is in excess of the 
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Wgvelertgth [A] 

Figure 1. Upper panel: Spectra for Channels 3 and 4 for all 1822 spectra analyzed in this 
work, stacked in order of Tqqes increasing upwards, with the scale shown on the left. Lower 
panel: Five RESIK Channels 3 and 4 spectra averaged over 2-MK intervals, with the averaged 
temperature Tqoes shown in each case (e.g. "5 MK" indicates spectra with Tqoes between 
4 MK and 6 MK). The Si XIV Ly-/3 line at 5.217 A has increasing intensity relative to the 
Si XIII Is^ — ls3p line at 5.688 A for increasing TgoeSi while the Si XIII line predominates in 
Channel 4 spectra at lower temperatures. Note that there is a short gap between Channel 3 
and Channel 4 spectra (4.86-5.00 A for on-axis flares). 



maximum of the theory curves. An explanation for this departure is discussed 
in Section [5] 

Shorter-wavelength, weaker lines of Si xiv are available for analysis in Chan- 
nel 3, viz. the Ly-S (Is — 5p) and Ly-e (Is — 6p) hues, but we defer this to a later 
work when the Si xiii atomic data become available. 
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Figure 2. Left panel: Si XIV Ly-/3 line flux [photon cm~^ s~^, estimated from the photon 
counts in the range 5.185-5.241 A with background in a neighbouring range subtracted] 
divided by EMqoes for spectra taken during the 21 flares listed in Table [l] plotted against 
Tgoes- 1669 spectra with more than ten photon counts in the Si XIV Ly-/3 line are included. 
Theoretical G(T) functions are plotted for photospheric (dashed-blue curve) and coronal (con- 
tinuous-red curve) abundances. Right panel: Histogram of values of j4(Si) derived from the 
Si XIV Ly-/3 line observations plotted in the left panel. The best-fit Gaussian curve indicates 
a mean abundance of j4(Si) = 7.92 ± 0.21. (The dashed-blue line indicates photospheric Si 
abundance, the continuous-red line the coronal Si abundance.) 



G(T) for Si XIII w3 band 5.66-5.71 A 

5 l ' 



Ag; w3 in Cli4 liistogram 




,[MK] 



7.5 8.0 8.5 
Log Abund [Ah=1 2] 



9.0 



Figure 3. Left panel: Si XIII Is^ — ls3p (w3) line flux [photon cm~^ s~^, from photon counts 
in the range 5.66-5.71 A with background in a neighbouring range subtracted] divided by 
BMqoes foi* the spectra taken during the 21 flares listed in Table[TJ plotted against Tqoes- 
1684 spectra with more than 50 photon counts in the Si XIII ui3 line are included. Right panel: 
Histogram of values of j4(Si) from the Si XIII ui3 line observations plotted in the left panel. The 
best-fit Gaussian curve indicates a mean abundance of j4(Si) = 7.89 it 0.13. (The dashed-blue 
line indicates photospheric Si abundance, the continuous-red line the coronal Si abundance.) 
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4. Flare Abundance of Si 

The analysis leading to the abundance of silicon in flare plasmas follows our 



procedure in previous work, see, e.g., Sylwester et al. (j2012p for sulphur lines in 



RESIK spectra. In brief, for a particular (ith) RESIK spectrum, a factor /i(Si) 
is estimated given by 

p. 

^'^^'^ ^ G(T,)EM, 

where Fi is the flux of the Si xiv Ly-/3 or Si xiii w3 line, the GOES temperature 
and emission measure at the time of the spectrum are Ti and EM^, and the contri- 
bution function [G] is evaluated at Ti and with an assumed Si abundance that is 
either photospheric ( |Asplund et al, 2009D or coronal ( Feldman and Laming, 200"0l ).| 



The measured Si abundance is then given by fi times the assumed abundance. 
A histogram is then constructed to show the distribution of the derived abun- 
dance values. On a logarithmic abundance scale, the distribution of estimated 
abundances has been found from our previous work to be close to Gaussian. The 
peak of the Gaussian distribution gives the most probable abundance estimate 
and the width the estimated uncertainty. 

This has been done for the Si abundances from the Si xiv Ly-/3 and Si xiii 
uiS lines, as is shown in the right-hand panels of Figures [2] and |31 We derive 
A(Si) = 7.93 ± 0.21 and 7.89 ± 0.13 respectively. The smaller uncertainty of 
the measurements of Si Xiii line, which is more intense than the Si xiv line for 
most of the temperature range, reflects the marginally tighter distribution of 
the observational points in Figure [H The slight departure of points at higher 
temperatures, already noted, does not have an appreciable effect on the Gaus- 
sian shape, in particular the peak and width of the curve. This is most likely 
connected with the accuracy of the atomic data in evaluating G{Te): the rate 
coefficients for excitation to the Si XIII n = 2 levels are from the distorted wave 



calculations of Zhang and Sampson (|1987p . but for higher n the rate coefficients 
are based on hydrogenic values. There are at present plans to incorporate into 
the CHIANTI data files much more precise rate coefficients from the i?-matrix 



calculations of Aggarwal and Keenan (j2010p . This will enable analysis of the 
other members of the Si xiii lines (Is^ — Isip, Is^ — ls5p) to be analyzed. There 
is, then, close agreement of the abundance estimates for these two lines from 
different ionization stages of Si. 



5. Discussion and Conclusions 

Figures [2] and |3] illustrate clearly that the flare abundance of Si from both the 
Si XIV Ly-/3 line and Si xiii w3 line is significantly higher than the photospheric 
abundance. With the latter taken to be A{Si) = 7.51±0.03 ( |Asplund a/., 2009[ ),| 
we find the enhancement of the flare value over the photospheric (called the FIP 
bias bv IFeldman et all [MIS)) to be 2.6 ± 1.3 (Si xiv) and 2.4 ± 0.7 (Si Xlll). 
Our Si abundance estimates are remarkably constant from flare to flare and over 
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the course of the development of each flare as shown by the smaU scatter of 
points in Figures [5] and [31 Table [T] shows that for the 21 flares analyzed the 
GOES peak emission had a range of a factor of > 400. This was found to be the 
case from RESIK measurements for the abundances of at least Ar and S (see 
[Sylwester et al. I2012p and from RHESSI measurements for the abundance of Fe 



dPhillips and Dennis, 20T2| ) 



Our measurements of flare Si abundances may be compared with coronal 
abundance estimates from other work, although such estimates are relatively 
few from flare spectra. Thus. FVeck and Parkinson] (| 1981 p used flare spectra from 
a crystal spectrometer on OSO-8 to give ^(Si) = 7.7g §, just consistent with 
our values from the Si xiii and Si xiv lines. Our values are higher than the 
two values (7.62, 7.52) measured from Si xiii 6.65 A line emission during flares 
seen with the Flat Crystal Spectrometer (FCS) on the Solar Maximum Mission 
determined bv lFludra and Schmelzl (|1995p . although it should be noted that the 
6.65 A line falls very near an anomaly in the FCS crystal reflectivity which may 
affect its intensity. 

The analysis of RESIK spectra during a low-activity period in January and 
February 2003 by |Sylwester, Sylwester, and Phiflips] dMni) was done by various 
methods including the derivation of the differential emission measure (DEM). 
The spectra chosen were from time intervals during the first quarter of 2003, 
grouped into five GOES level ranges, from A9 to B5. The derived DEM, from the 
Withbroe-Sylwester method ( [Sylwester, Schrijver, and Mewe, 1980[ ), generally 
had a bimodal distribution, with a low-temperature component (2-3 MK) and 
with a high-tcmperature component (6-9 MK) which becomes more significant 
with rising GOES level. An inverse relation of the temperature of the hotter 
component with GOES level was found, but with A(Si) almost constant, at 
« 7.4, about a third of the value obtained here. 

Spectra from the Hinode Extreme-ultraviolet Imaging Spectrometer (EIS) have 
been analyzed to give Si/S abundance ratios from quiet-Sun and active region 
spectra - the stigmatic nature of the spectrometer allows individual coronal 
features to be examined. Taking S to be a high-FIP element, i.e. one with 
FIP > 10 eV, the Si/S abundance ratio therefore gives the FIP bias of Si, 
defined by [Si/S]cor/[Si/S]phot- The recent analyses bv lFeldman et a/.l(|2009p and 
IBrooks and WarrenI (|201ip have used the fiux ratio of Si x and S x lines in EIS 
spectra to find the Si/S abundance ratio and so the FIP bias of Si. Abundance 
studies using extreme ultraviolet spectra have a number of difficulties, which are 
enumerated by IFeldman et al.\ (j2009p . Among the most significant is the very 
crowded nature of this spectral region compared with, say, the RESIK soft X- 
ray range, so that the possibility of line blends is often high. With this proviso 
in mind, IFeldman et al.\ (j2009p obtain values of Si FIP bias that range from 2.3 
(quiet-Sun limb) to 5 (quiet-Sun limb but with "moderate activity"), and with 
active region values of 4.6-4.8. IBrooks and WarrenI (|201ip observe the active 
region outflows which have been widely observed with EIS spatially resolved 
spectra, finding the FIP bias to be 2.5 to 4.1. FIP bias values for the slow solar 
wind measured by these same authors using data from the SWIGS instrument 
on Advanced Composition Explorer (AGE) are found to agree with the latter 
EIS result for a particular active region apparently magnetically connected to 
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the location of ACE, so indicating the connection of the slow solar wind with 
the active- region outflow observed by EIS. 

The RESIK observations of solar X-ray emission during a period of low activ- 
ity ( [Sylwester, Sylwester, and Phillips, 2010[ ) include both S and Si lines from 
which the Si/S abundance ratio can be derived. Although the Si abundance was 
nearly constant for the five GOES ranges, the S abundance varied from A{S) = 
6.75 to 7.25, inversely related to GOES activity level but positively related to 
the temperature of the hotter component of the bimodal DEM distribution. 
The Si/S ratio increases from 1.6 to nearly 5 for increasing temperature of the 
hotter component, similar to the range of Si/S ratios measured bv lFeldman et aH 
(|2009p . For comparison, the flare Si abundance from the data discussed here and 
the flare S abundance from [Sylwester et a/.|(|20f^ give a Si/S abundance ratio 
of 5.43 3 taking into account the uncertainties of both the Si and S abundance 
measurements. The Si/S ratios of [Sylwester, Sylwester, and Phillips] (|2010p may 
fit a picture in which the Si/S abundance ratio steadily increases in active regions 
with increasing activity to fully fledged flares as studied here. 

In summary, the present analysis of Si xiii and Si xiv lines seen in RESIK flare 
spectra give estimates of the Si abundance of A{Si) = 7.93 ± 0.21 (Si xiv Ly-/3 
line) and ^(Si) = 7.89 ± 0.13 (Si xiii wS line), with very small variation from 
flare to flare or times within flares, in approximate agreement with measurements 
by IVeck and ParkinsonI (|198ip . These abundances are a factor 2.6 ± 1.3 and 
2.4 ± 0.7 greater than the photospheric abundance ( [Asplund et ai, 2009 ). The 



Si/S abundance ratio for flares as observed by RESIK is large compared with 
values for non-flaring active regions, suggesting that the ratio depends on the 
degree of activity. 

In future work, we will use improved atomic data, when available in the 
CHIANTI package, to obtain Si abundance estimates from Si xiii lines seen in 
RESIK Channel 3 spectra. 
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